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Abstract 

Using exact results from the theory of completely integrable systems of the 
Painleve/Toda type, we examine the consequences for the theory of polyelec- 
trolytes in the (nonlinear) Poisson-Boltzmann approximation. 

1. Introduction 

A polyelectrolyte is a macromolecule with a large number of ionizable groups 
which in a solvent becomes highly charged |1|, 0]. This highly charged structure 
is referred to as the polyion. If the polyion has charge Q = Ze, then in solution 
there are Z counterions. (We take the polyion to have negative charge and thus the 
counterions have positive charge.) Typically the solvent also contains a salt and so 
the total number of counterions is Z plus the number of positive ions from the salt. 
The negative salt ions are called colons. Typically the excess salt is a 1-1 salt, e.g. 
NaCl, though other salts are considered. (MgCl2 is an example of a 2-1 salt.) 

One class of polyelectrolytes that is widely studied are those consisting of slender, 
rod-like particles. For example, the tobacco mosaic virus, a polyelectrolyte, has a di- 
ameter of about 18 nm with a length of 300 nm. For such systems the idealized model 
is an infinite cylinder of radius a of uniform linear charge density with counterions 
and colons treated as point particles. This neglects several important effects such 
as the interaction between polyions, the flexibility degrees of freedom of the polyion, 
finite size effects, to name but a few. Nevertheless, this model has been extensively 
studied (for reviews see Q). In this idealized model there are two additional 
length scales: the electrostatic length scale (Bjerrum length) is = e'^/ekBT {e is the 
solvent permittivity) and the inverse Debye-Hiickel screening parameter = Stt^b/ 

^Dedicated to Benjamin Widom on the occasion of his seventieth birthday. 



where I = ^jflj is the ionic strength with rtj the molar concentration of ion j of 
(integer-valued) charge qj with the convention that qj < for counterions and qj > 
for colons. From these two lengths we form the dimensionless Manning parameter 
B H ^ = ^-b/^ where b is the average spacing between charges on the polyion and the 
dimensionless distance parameter r = kR where R is the cylindrical distance variable. 

A mean field theory approach to this model results in the Poisson-Boltzmann 
(PB) equation; that is, the electrostatic potential \& is assumed to satisfy the Poisson 
equation of electrostatics with the density of various ions being given in terms of 
Boltzmann factors |P, 0, H, ^. Introducing the reduced potential y = /ksT the 



PB equation becomes 



-liV 



In terms of r and for cylindrical geometry the PB equation is, explicitly for the cases 
of 1-1 and 2-1 "excess of salt" 0], 



cPy Idy _ { sinhy 1-1 salt, , 



dr"^ rdr M (e'^ - e"^) 2-1 salt. 



One boundary condition for ( p..l| ) is obtained by applying Gauss's Law at the 
surface of the polyion 

lim = -2i. (1.2) 

(We use the same symbol a for both the polyion radius and the dimensionless polyion 
radius.) The second boundary condition depends upon whether the system is closed 
or open. For a closed system one requires that the electric field vanish at some finite 
distance whereas for an open system with a finite concentration of added salt one 
requires that 

y{r) — > as r CX3. (1.3) 

We will consider only the case of an open system. A further simplification, and one 
we will also assume, is to take the polyion radius a —>■ which means we model the 
polyion as a line charge. The mathematical problems that we address in this paper 
are now well formulated: To solve (|1.1|) subject to the boundary conditions (|1.2|) and 

(ra. 



2. Exact Solutions for a 1-1 and 2-1 Salt 



First note that the solutions to the linearized versions of ( |1 . 1| ) satisfying ( |1.2| ) and 
3D are in both cases the familiar Debye-Hiickel solution 

VDHir) = 2^Ko{r) (2.1) 
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where Kq is the modified Bessel function. The exact solution for the 1-1 sah is 



m 11, 12 



2/11 (r) = 21ogdet(/ + A/s:) -21ogdet(/- Air) 

oo \2j+l 

= -Tt(k''+' 



,=0 2J + 1 

where K is the integral operator on with kernel 

exp (— r/2(x + 1/a;)) 



(2.2) 



and 



X + y 



A = — sin — for E < I. 

TT 2 ^ - 



Note that Tt{K) = Ko{r) and that Tt{K^] 
been proved that ||10|, |l3l 



2/11 (r) = —2^ log r + 6C, log 2 + 2 log 

where F is the gamma function. 

The exact solution for the 2-1 salt is 



0(e as r — i> oo. For r — > it has 

r(¥) 



F(i^ 



o(l) for ^ < 1 



(2.3) 



2/21 J 



log det (/ - XKs) - log det (/ - Ai^s) 



A^' 



E i Tr (K 



j=0 



Tr (iC^ 



(2.4) 



where Kj (j = 2, 3) are integral operators on TV~ with kernel 

uj^-^{l - u) exp [-(r/2V3) ((1 - uj)x + (1 - cj-i)^-^)] {-cux + y)~^ 
+ cu2{^^-2)(l - uj^) exp [-{r/2VS) ((1 - uj^)x + (1 - uj-^)x-')] {-uo^x + y) 



-1 



for i< -. 



Here u = e^'^*/^ and 

A = 2sm(fK+l/4)J-l. A.= ^ 
Again note that 

TT{K2)-TTiKs)=6Koir). (2.5) 

An elementary calculation produces a simpler representation than the definition for 
the trace of the square of the operators: 



exp (-r/2(xi + 1/xi + X2 + I/X2)) 



Xi + 0:10:2 + ^2 



dxidx2 . (2.6) 
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The computations ( |2.5| ) and ( |2.6| ) lead one to suspect that for all postive integers 
n the quantities Tr(fr2 ) ~ Tr(i^'3), as functions of r, are monotonically decreasing. 
Indeed, one can derive an alternative matrix kernel representation of the operators 
Kj such that the monotonic decay is manifestly clear. 
For r — > it has been proved that O, IP 



2/21 (r) = -2^ log r + (2 log 2 + 3 log 3) ^ + log + o(l) for ^ < 1/2. 

^ V 3 j ^ V 3 j 

(2.7) 

Higher order terms in both ( |2.3| ) and ( |2.7| ) can be computed from use of the differential 
equations and all additional constants appearing can be expressed in terms of the 
quantities above. 

3. Asymptotics at the Critical Value of ^ — Counterion Condensation 

The Oosawa- Manning arguments fl], ^, ^ for counterion condensation are well 
known and need not be repeated here. These arguments, when applied to a 1-1 salt 
and 2-1 salt, predict critical values of ^ at 1 and 1/2, respectively. Recall that in 
the theory of counterion condensation, the meaning of ^ > is that the average 
charge spacing h on the polyion is increased by counterion condensation (onto or 
around the polyion) until ^ = is achieved. Thus one must distinguish between the 
stoichiometric value of ^ computed using only the charge groups on the polyion and 
the lower value of ^ achieved through counterion condensation. 

Mathematically, the critical value of ^ is seen through the qualitative change in 
the small r asymptotics (p.3|) and (|2.7|). For the 1-1 salt case we have as r ^ 



m 13 
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exp (-2/11 (r)/2) = -Ul, - ^ {snl - Snj + 4fii - l) + Oir'nl) (3.1) 

where 

ni = ^]l(r) =log(r/8) + 7 

and 7 is Euler's constant. Thus the potential yu{r) develops an additional log log r 
singularity at the critical Manning parameter. Similarly for the 2-1 case as r — * 

EH 

' 2„ 2^ 



exp (-2/21 (r)) = + ^ (^2 - 3^2 + g) + 0(r^f^^) (3.2) 



where 



^2 = ^^2 (r) = log(r/8) + 7 + ^ log 2 + ^ log 3. 



For ^ > the above solutions are no longer physically valid. For example, the 
Boltzmann factor exp (— ?/ii(r)) will become negative for small enough r. 
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4. Exact Electrostatic Free Energy for 1—1 and 2-1 Salt Cases 



At constant temperature and pressure the free energy is the work done in placing 
charges on the polyion. This is the familiar "charging process" and one imagines an 
increment of charge dq placed at the surface of the polyion so that the infinitesimal 



work done is ^{a)dq where the electrostatic potential is evaluated at R = a [[15], [16 
Thus 

„1 rQ 
w^'- = ^(a) dq 
Jo 

is the free energy associated with a single line charge. Our solutions yn and ?/2i are 
of the form 

?/(r) = -2^1ogr + ?/o(0 +o(l) as r ^ 
and are for the limiting case of a 0. We take, therefore, for the value of the 
e^{a) /ksT the quantity — 2,^ log(Ka) + 2/o(0- Re-writing the above expression for 
w^^ in terms of dimensionless quantities and multiplying the result by Np^ the total 

number of polyions in solution of volume V ^ we obtain for the free energy, W^^^ of 
the entire solution fH, Bl 



lyel 



VkBT 



Up Mog(Ka)^ - ^ J^yoiO di' 



where Up = Np/V is the polyion concentration. Using the expressions (|2.3| ) and (|2] 
we calculate /(^) in these two cases: 

f^jn = (log(/s:a) -31og2)e + - / log ) ^ ' ,[ d^ 
jin^j V 6V ; 6 ;s ^ ^r(l/2 + f/2) ^ 

= (log(/ta) - log2 + 7) e - V 

V &v ; & ^^22"-i(2n + 2)!^ 



(log(fi:a) + l-31og2) + 21og 



2 



2^ r(¥)r(¥) , 4^ g(¥)g(¥) uu 

X A . X 1 o 3 1 oA . 1 /"^ r (2/3 - ^73) r (1/3 -2^73) 

/2i(^) = log(/ta) -log2- -log3K + - / log— r^^ ^^r^ — t^, 

&v ; & 2 ^ ^Jo ^ r(l/3 + f/3)r(2/3 + 2f/3) ^ 

= (log(«:a) - log2 + 7) e + ^ (^^'^1/3) - ^^^^2/3)) e 

(V^^'Hl/3) + ^(^^2/3)) e + 0(a. (4.2) 
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Free energies 




Figure 1: The free energies fu and /21 as functions of the Manning parameter ^ < 
with K,a = 1. At the critical value, /ii(l) = logna + 0.674563 and /2i(l/2) = 
i log Ka + 0.393915. 



In the above ijj{x) = r'(a;)/r(a;) and ijj^"'\x) is the n derivative of ijj. The function 
G{s) is the Barnes G-function [|17| which is an entire function of s and is defined by 

G{s + 1) = (27r)^/2 exp {-s/2 - (1 + 7).V2) 11 (l + ^) ^^P ['^ + ^ 

The G-function satisfies the functional equation G{s + 1) = r(s)G(s) and has the 
special value G'(l) = 1. It arises in the present context through the integral 

r log r(a; + l)dx = ^ log 2n - lz{z + 1) + z log T{z + 1) - log ^(z + 1). 
Jo 2 2 

(Clearly one can also express /21 in terms of the G-function.) 

In the small ^ expansion for both fu and /21 the linear term is the contribution 
from the Debye-Hiickel theory. Observe that since the volume dependence resides 
solely in the logna term, such derived quantities as the osmotic pressure (which is 
expressible in terms of a first derivative of W^^ with respect to V) will be identical 
to that derived from the Debye-Hiickel theory. The mean activity coefficients 7±, 
expressible in terms of the function /, will have different numerical values from the 
Debye-Hiickel theory, however their dependence on Up will be identical. Though this 
has been understood and used by previous workers, we stress these facts since it 
emphasizes the wider validity of the linear theory than one might a priori expect. 
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Both fii and /21 are singular at the critical Manning parameter and have leading 
singularities of the form {C,c ~ log(.^c ~ 0- The graphs of fu and /21 are shown in 
Fig. 1 where for convenience we have set na = 1. 

5. Partial Equilibrium Structure Factors 

Within the Poisson-Boltzmann approximation the density distributions of coun- 
terions (+) and colons (— ) are given by 

n±(r) = n± exp {^qjy{r)) 

where n± are the concentrations at infinity. Light scatters from the local concentration 
fluctuations of the counterions and colons and the observed scattering intensity in 
the static approximation is expressible in terms of the partial structure factors which 
are, in the cylindrical PB approximation, the (two-dimensional) Fourier transform of 
{n±{r)/n± — 1) [l^, Thus we examine, for the 1-1 salt. 



S±{q, = -^(e^^ - 1) = 27r Jo{qr) (e^^^") - l) r dr (5.1) 

where q is the dimensionless wave number and Jq is the Bessel function of zeroth 
order. For the 2-1 salt we replace — »• e^^ in the above expression. 
At g = 0, 

c (nn _ / 47re(l + e/2) 1-1 salt, 

^+^^'^> ~ \ 87re(l + e/2) 2-1 salt, ^^■^> 
c _ f-4<(l-^/2) 1-1 salt, 

•^-^^'^^ ~ 1-47re(l-0 2-1 salt. ^^'^^ 



The differences S+{0,^) — S'_(0, ^) follow directly from the differential equations ( |1.1| ) 
and the boundary condition (|1.2| ). The individual terms require additional identities 
(for 1-1 case see eq. (6) in [|lll). In the Debye-Hiickel approximation ( p^TD , S±{0,^) 
is ±4^ for the 1-1 salt and similarly for the 2-1 salt. For the 1-1 salt, S'_(g, ^) = 
S+iq, -0- 

To understand further how these results differ from the Debye-Hiickel approxima- 
tion, we expand the exponentials in ( |5.1|) and use the expansions (p^) and ( |2.4|) to 
deduce the expansions: 

00 

For the 1-1 salt we find, for example. 



1 + g2 



7 



1-1 Salt 
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Figure 2: The normalized moments /i± as functions of the Manning parameter ^. 



log g/2 + J(g/2)2 + l 



(g/2)^(g/2)2 + l 



32n 



{Xi +X2 + X3)(l/Xi + 1/X2 + 1/0:3) 



^0 JO (xi + X2){X2 + X3){X3 + Xi) 

(Xi + 1/xi + X2 + 1/X2 + X3 + I/X3) 



[4g2 + (xi + 1/xi +X2 + 1/X2 +X-i + I/X3)] 



3/2 



dxidx2dx3 



From ( |5.2D it follows that S'-|-^3(0) = ±47r'^/3. In general S±^2p+i and 5'-|-_2p+2 will 
involve the Tt{K'^^'^^) {i = l,...,p) and they will have, in the complex g-plane, 
branch point singularities at g = ±ip. Since A = O(^), ^ — >• 0, we see, for ^ -C 1, that 
the Debye-Hiickel approximation of taking the first term in the sum ( [5.4| ) is valid for 
bounded q. Similar expansions can be derived for the 2-1 salt where we note that 
S±^i{q) will again be the Debye-Hiickel term (see ( |2.5| )) but now S±^2{q) will contain 
irreducible "two-particle" effects (see (p.6|)). 

For any scattering intensity I{q), one experimentally accessible quantity is the 
normalized second moment n in the small q expansion 



m 
m 



l + /ig2 + 0(g' 



For the partial structure factors S± one similarly defines /i±(0 and in Fig. 2 we 
graph fi±{^). We have defined the inverse correlation length by the distance to the 
nearest pole in S{q). (The choice of the dimensionless wave number q = k/n fixes 
the location to For the Debye-Hiickel approximation /i± = 1, in contrast with 
the PB equation where (yU_) decreases (increases) with increasing ^. The partial 
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structure functions themselves are shown in Figs. 3 and 4 for various values of ^. Note 
that with increasing ^ there is a significant increase in 5'+(g, ^)/S'+(0, ^) for large q 
whereas S'_(g, ^)/S'_(0,^) shows only a slow decrease for increasing ^. 

To determine the large q asymptotics of S±{q,C,), we use the fact that if 

/(r) ~ as r — > 0, 

then 

with C(A) = -2^+^ sin(7rA/2) 1(1 + X/2f. From Q and the above fact we deduce 
for the 1-1 salt that as g ^ cxd, for fixed ^ < 1, 

S±i<l,0-^ (5.5) 

where C± = exp (±?/o(^)) C(=f20- At the critical value ^ = 1 the short distance 
asymptotics of e^" are no longer of the above form (see ( |3.1| )). For this case we find 
as q —>■ oo 

SMSJO,!) = 1^-^ 1 +0( ' 



3 logg 3 log q log g 
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S+(q)/S+(0), 2-1 Salt 



0. 6 
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S- (q) /S- (0) , 2-1 Salt 



0. 4 




Figure 4: The structure factors S±{q,^) / S±{Q,^) as functions of q for the 2-1 salt. 
Also plotted is the Debye-Hiickel approximation. 



5_(g,l)/5_(0,l) 



21ogg ^ 2(2 log 2 



+ lower order terms. 



Similar large q expansions hold for the 2-1 salt. 

In critical light or neutron scattering from simple fluid or magnetic systems, the 
large q behaviour of S{q) defines the critical exponent r] For mean field 

theories of critical scattering, 77 = 0. We remark that even though the PB equation is 
a mean field theory, we have a nonzero "r/" [t] = ±2^, ^ < 1, for the 1-1 salt) which is 
a reflection of the fact that the short distance potential is the bare Coulomb potential. 
As the critical value of the Manning parameter is approached, the Poisson-Boltzmann 
theory predicts an enhanced scattering at large wave numbers from the concentration 
fluctuations of the counterions while at the same time predicting little change in the 
scattering from the concentration fluctuations of the colons. As is the case in critical 
scattering 24 1, the measurement of these effects may well prove difficult. 
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